Resonant Optical Nonlinearity of Conjugated Polymers 
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When the energy of a pump wave is in resonance with the exciton creation energy, the electric 
susceptibility of a conjugated polymer in response to the probe wave is altered by the exciton gas. In 
this paper, we calculate the dependence of this change on the the exciton populations by the equation 
of motion (EOM) method. The magnitude of optical nonlinearity is also influenced by ambient 
temperature, by the extent of exciton wave functions, and by the strength of electron-electron 
interaction. All of these factors can be easily incorporated in the EOM approach systematically. 
Using the material parameters for polydiacetylene (PDA), the optical Kerr coefficient 712 obtained 
is about 10~^ cm^/W, which is close to experimental value, and is four orders of magnitude larger 
than the value in nonresonant pump experiments. 



PACS numbers; 78.66.Qn 



I. INTRODUCTION 



As a class of materials promising for applications in 
all-optical devices for communication and data process- 
ing, conjugated nplymers have been a subject of great 
research interest.El Conjugated polymer such as PDA 
has long been recognized to exhibit large nonresonant 
third-order optical nonlinearitiesfl Most of the theoret- 
ical works on the nonlinear optics of these materials 
are based on the perturbative sum-over-state formulaEl 
that suits nonresonant pumping well (but is inappropri- 
ate for resonant pumping). A fair agreement between 
experiments and theories has been achieved for a vari- 
ety of nonresonant spectra, including two-photon absor 
tion, third-harmonic generation, and electroabsorption 
In addition, EOM approach has also been employed in 
these studiesJj The EOM for polarization, which is sim- 
ilar to the semiconductor Blach equation originally used 
for inorganic semiconductors^, is a general tgol whose va- 
lidity goes beyond the perturbative regimeQ and will be 
used in this paper. 

For conjugated polymers under resonant pumping, 
most of the optically excited electrons and holes remain 
bound at room temperature and form an excitorijgas, in- 
stead of a plasma as in the IH-V semiconductors. 3 This is 
because the magnitude of the exciton binding energy in 
conjugated polymers (of the order of 1 eV) is much larger 
than that in inorganic semiconductors (a few meV). This 
presents a difficult situation for the calculation of reso- 
nant nonlinearity. Even though there are a lot of exper- 
imental works, few theory has been devoted to this sub- 
ject. A simple and heuristic explanation of this nonlinear- 
ity based on phase-space-filling (PSF) has been proposed. 



and the result agrees quantitatively with experiments. □ 
The esssence of PSF is that when there is a finite concen- 
tration of excitons, the phase space for further excitons 
to form is reduced because of the exclusion principle, and 
henceforth the probe absorption signal is reduced. The 
limitation of PSF is that it cannot predict the whole spec- 
tral response of the electric susceptibility, nor the effects 
of temperature, strength of electron-electron interaction 
... etc. 

Motivated by the desire to have a more accurate theo- 
retical tool to understand this phenomenon, we present a 
detailed analysis of the nonlinear susceptibility based on 
EOM. We will utilize the technique developed by Haug, 
Koch, and Schmitt-Rink to make a connection between 
exciton populaticms and the electron populations that ap- 
pear in the E0MJ13 However, since the link is suitable for 
dilute exciton gas only, we shall concentrate our study 
in this regime. Furthermore, from the comparison be- 
tween the lineshapes of timc-rcsolved and cw photolumi- 
nescence spectra, the typical intraband scattering tirne. 
is estimated to be less than picosecond in polymers. Ell 
Therefore, it is possible to consider a simpler situation 
where the exciton gas is in quasi-thermal equilibrium, 
whose population is determined by the intensity of the 
pump wave. Because of these approximations, phenom- 
ena such as exciton-exciton interaction, excited-state ab- 
sorption, and off-equilibrium momentum distribution are 
not considered. 

In this paper we have done extensive studies on the 
optical nonlinearity of PDA. There are several advan- 
tages in choosing this material: First, it has the largest 
nonresonant third-order nonlinearity of all polymers and 
very large resonant nonlinearity. Second, it is one of a 
few conjugated polymers that can form high quality sin- 
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gle crystals and is more amenable to theoretical analy- 
sis. Third, its chain-to-chain distance is large because of 
large side groups and therefore the inter-chain interac- 
tion is less important. Various aspects of the nonlinear 
spectra are investigated. By choosing a reasonable value 
of the Coulomb interaction strength, the resonant optical 
Kerr coefficient 712 being calculated agrees very well with 
experiments. Furthermore, the effect of electron-electron 
interaction on the height and position of the exciton ab- 
sorption peak is studied. We also study the influence 
of temperature, as well as the relative populations of sin- 
glet and triplet excitons, on the probe spectra. This work 
provides a systematic analysis of the influence of various 
microscopic parameters on the optical nonlinearity and 
may serve as a guide for the search of optical materials 
with larger resonant nonlinearities. 

This paper is organized as follows. In section 11, the 
equation of motion for polarization in conjugated poly- 
mers is derived. In section III, the Coulomb potential ma- 
trix elements and electron occupation numbers are cal- 
culated. The numerical results are presented in section 
VI, and Sec. V is the conclusion. 



II. EQUATION OF MOTION FOR 
POLARIZATION 

A conjugated polymer is a macro-molecule with many 
electronic and ionic degrees of freedom. There have been 
several attempts to include both types of degrees of free- 
dom in calculations on the electronic and optical prop- 
erties of conjugated polymers. This is a very difficult 
task. To date, an exact solution for chains longer than 20 
sites is still lacking.t3 It is especially challenging to calcu- 
late resonant optical nonlinearity by including the effects 
of electrons (including excitons) and phonons simultane- 
ously. In this paper, we choose a modest approach and 
focus only on the electronic contribution to the nonlinear- 
ity. By doing so, we will not, for example, be able to pro- 
duce satellite phonon peaks next to the main resonant ab- 



sorption peakJ13 Our main goal is to calculate the optical 
nonlinearity due to the exciton gas, which is derived from 
the dependence of the magnitude of the main absorption, 
peak on the exciton populations. In fact, Greene et auu 
have demonstrated that the PSF model, whose origin is 
purely electronic, can explain major features of the spec- 
tra well. Furthermore, we consider only the dynamics of 
TT-electrons; the c-electrons are tightly bound to the ions 
and has little influence on the dynamical response. How- 
ever, they do contribute to the renormalization of the 
interaction between 7r-electrons, and between ions and 
TT-electrons. These effects appear implicitly through the 
parameters in the 7r-electron Hamiltonian. 

Our calculation is based on the SSH-like Hamiltonian 
with electron-electron interactions, H = Hq + Hi + H2, 
where 

-ff = ^ {^ckdlka^cka + ivkal^^ayk^^ , (2.1) 
ka 

TT \ \ ■\)k\k2k'^k4^ t t 

^1 ~ '^AiA2A3A4"Aifcicr"A2A:2cr'"A3fc3cr'"'A4fc4CT, 
A1A2A3A4 kik2k^k4^ 

(2.2) 

and 

H2 = -E{t) ^ (^tJ-vcik)al^^acka + lJLcv{k)a\j^^avka^ ■ 

ka 

(2.3) 

In these equations, ec,vk are the energies for the 
dimerization-induced conduction and valence band, 

yfeifc2fc3fc4 g^j,g ^jjg Coulomb potential matrix elements 
A1A2A3A4 ^ 

(A = c, u), and /iAA'(fc) are the dipole matrix elements. 

H2 describes the coupling between polarization and the 

prohe field E{t) along the chain. The influence of the 

pump field will be accounted for when calculating the 

conduction electron populations in the next section. The 

potential matrix elements are 

(2.4) 



V 



ki k2k3k4 
Ai A2 A3 A, 



:= / Mr'n,kMn2k2{r')Vee{r-r')^x,kArl^x,k. 



where Vee is the screened Coulomb interaction between 
the TT-electrons, and '^\k are the Bloch states solved from 
Hq. The usual practice is to neglect, from the very begin- 
ning, the terms that do not conserve electron numbers in 
each band, which include half of the sixteen terms that 
do not have equal numbecs of c indices and v indices, plus 
V^icc ''''''' and V^i^^'^'^'^'^H This is reasonable in metals or 
inorganic semiconductors where the Coulomb interaction 
is weak, but is not necessarily valid in conjugated poly- 
mers where the Coulomb interaction is much stronger, as 
indicated by the large exciton binding energies. There- 
fore, these terms will be kept in the derivation until they 



are proven negligible. It will be shown later that the 
terms without equal numbers of c and v indices indeed 
make no contribution, but the V^l'^f^''^ and V^i,^^^^'"' 
terms cannot be ignored. In fact, unlike all the other 
terms, these two terms do not conserve the electron and 
hole spin individually, thus are essential to the lifting of 
the four-fold degeneracy in the spin subspace of the ex- 
citon states. 

In the following, we derive the EOM for pka = 
{a\j,^ar,kiy) ■ The total polarization (P), which is equal 
to J^ka(t^'"cik)pka- + fJ-cvik)pl^), Can be easily obtained 
by integration over all fc's. In general, terms of the 



2 



form Efc,fe,fe3i:.'"^A^ \l\l'\^\ikia^\2k2<T''^>'aka<T>axkcr are mean values of the product of four operators can be fac- 
encountered in the derivation. By using the RPA, the torized, 



(«AifeiO-«Lfe2<T'"^3fc3'T'aA4fc4<T> (^li fei <T« ^4 fcl <T ) fea <t' « ^3 fea 't' ) ^fel fe4 '5*2 fes 

~ (»AifciCr"A3fcl<T)(al2fc2Cr"A4fe2<T)4lfe3 42fe4'5o 

After a straightforward but tedious calculation, we have 

. dpka 



ih- 



dt 



fc' 
fe' 

+ E 



i^vkcT ^cka) 



{j^vka T^cka^^ 



(2.5) 



(2.6) 



*c vka'^c,vka) are the one electron occu- 



where Ucyka 

pation numbers in conduction or valence band, V 



kk' 

Ai A2 A3 A 

and V^'^^^^^)^ are abbreviations for the Coulomb terms. 



V 



kk kk 
Ai A2 A3 A ' 



and the exchange terms, VxixlxsX^ respectively. 
The terms that are quadratic in pka are neglected in 
Eq. ( |2.6| ), because only the linear response of the probe 



wave is considered. Also, relations such as V^^^^ — V^^^^ 

)kk' 

vcvc ■ 



)kk 



and Vi'^™ 



symmetry between conduction band and valence band in 

[(Ccfe - ei,fc ) - W - 17] pka {uj) 
= E{uj)^J.yc{nvka - ricka) 



the present model. Note that the potential matrix el- 
ements with unequal numbers of c and v indices result 
in the last line of Eq. (2.6). These matrix elements are 



multiplied by terms quadratic in electron or hole occu- 
pations, which under resonant pumping may be large. 
Notwithstanding, because of inversion symmetry, it can 
be shown that for thermalized exciton gas, these terms 
actually have no effect on the dynamics (see Appendix 
A). Consequently, Eq. (2^) becomes 



-E 

k' 

+ E 



fc' L 



'^vvvv ^cvcv^ {n-uk'a- 'T'cfc'cr) ^ ] (^vvvv '^cvvcj ('^ufc'cr' l^ck'a') 

a' 



Pka{uj) 
{j^vka ~ T^cka')-! 



(2.7) 



where pj^^ = {pka -\- P-ka) 1'^, and a damping term has 
been added. Notice that p\^{—uj) ^ pkai^) because pka 
is not Hermitian. Therefore, Eq. ( ^.7[ ) has to be solved 
in conjunction wit h th e equation satisfied by p^^, which 
is similar to Eq. ( ^.7[ ) but with uj replaced by and 
Pka replaced by p^^. The meaning of the various parts 
on the right hand side of the equation is explained be- 
low: The first square bracket, after being summed over 
/c', contributes to the band gap renormalization. The 
magnitude of renormalization depends on the strength of 
inter-electron interaction, as well as the electron popula- 
tions. Inside the second square bracket, the V^^™ term is 
most crucial to the formation of excitons; the V^^^^ term 
is related to the singlet-triplet splitting of the exciton lev- 
els and leads to an unusual coupling between positive and 
negative frequency components of the polarization. For 
the exchange terms, it will be shown in the next section 



that Vy^y^ = Vcvvc because of charge neutrality condi- 
tion. Therefore, there is no exchange effect in the first 
square bracket. We can also show that V^^uc ^^'^ '^vvcc 
are simply constants and can be treated as corrections 
to V^cct; and Vy^^^. The only unknowns in Eq. (2/7) are 
{pkai'-^)Tpkai~^)}- AH of the other quantities, including 
the potential matrix elements and the electron popula- 
tions, can be obtained given physical conditions such as 
the strength of electron interaction, the intensity of pump 
wave. . . etc. This is derived in the next section. 



III. POTENTIAL MATRIX ELEMENT AND 
ELECTRON POPULATION 
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A. Potential matrix element 



A natural choice for the interaction potential between 
electrons is Vb/|r — r'|, where Vq is given by e^/eao, e 
is the intra-chain dielectric constant, and oq is the av- 
erage distance between neighboring sites. The position 
r, being dimensionless now, is measured in units of cq. 



To calculate V 



kk' 

Ai A2 A3 A4 ' 



the unperturbed eigenstates are 



expanded by localized Wannier functions Jl^ 



/ M 

1 = 1,2 \^ m=l 



(3.1) 



where M is the number of unit cells. The 
total chain length is 2M. Defining = 

1 /2 

[e*'"'(to cos/c — i(5tsinfc)/ecfe] , where to — {—lySt is 
the hopping amplitude between neighboring sites, then 



V 



kk' 

Ai A2A3A4 



V2iq) 



U-ul(fc) Uv2{k) 



Ck a 

-Q Ck 



(3.2) 



When calculating the matrix elements of V"ee(r — r')j only 
the integrals involving Wannier functions at the sarpa 
site are kept (zero differential overlap approximation) .113 
To improve upon this, we need to know the shape of 
atomic orbitals, which will not be considered here. For 
the same site, there is a finite on-site energy Uq = 
1/2 / d^rd^r'W*{r)W*{r')Veeir - r')l^(r')l^(r). Since 
the exact form of the Wannier function is not known, 
Uq is treated as a parameter independent of Vq. So our 
choice of the potential is essentially of the Pariser-Parr- 
Pople form. Defining 



rn 

^2(«) = T?E^-(2m)e-2™^ (3.3) 



then a straightforward calculation gives 



"Aij (^)«A2j (^')wA3j {k)u\,j {k') 



-f V^{q) [e»«<^2(fc)<^i(fc')"A3l(fc)«A.2(fc') + e-'?<^l(fc)u*^2(fc')"A32(fc)^/Aa(fc')] , 



V 



kk' 



A1A2A3A4 - 2^2(0) 



Vm K2(fc)«A2l(fc')«A3l(fc')^A42(fc) +<,i(fc)<,2(fc')"A32(fc')«A4l(fc)] , 



(3.4) 



where q = k — k' . We require Vi(0) + ¥2(0) = be- 
cause of the charge neutrality condition. After com- 
bining Eq. (3.4) with Eq. (3.2), it can be shown that 
Vl^L - ViKc = (^2(0) + Fi(0))/2 = 0, and V^^:, = 
Vfe£^ = (1/2(0) - Vi(0))/2. We emphasize that this 
nonzero exchange correction is crucial to the reduction of 
the exciton binding energy when Uq is tuned to a larger 
value while keeping Vei fixed. The same behavior is ob- 
served in Abe's paperEj, w^re their concern is the energy 
spectrum of the excitons.llZl 



tion of an exciton at the n-th bound state with angular 
momentum {jm}. Helped by this relation, we can write 
n-cvk in terms of exciton operators as follows: By using 
an identity operator 

^ " W'^^ " '^ika'^^ka) °^"fc'^"L<T' (3-6) 

K<T k(7 



where N is the total number of electrons and N,, is the 
number of conduction electrons, we have 



B. Electron population 

The electron population can be linked to the exci- 
ton population by using the method developed by Haug, 
Koch, and Schmitt-Rink. Their relation is derived be- 
low. Firstly, the connection between electxon operators 
and an exciton creation operator ej^ ial^l 



''cki(Ji^vk20 




m\c7ia2) 



njmK ' 



(3.5) 



where {jm\aia2) is the Clebsh-Gordon coefficient, k = 
(fci + ^2)72, K = ki — k2, and 4>njm(k) is the wave func- 



^ A;2cricr2 



2MN, 



Y 'Y^*r^jrnik)(t)n'3rn{k)elj^j^en'jmK- (3.7) 



nn'jm k2 



For a dilute exciton gas, the exciton number 
{'^li^niK^n'jmK) — gnjmK^nrd^ where gnjmK IS the ther- 
mal distribution of excitons.li2l Therefore, 



= E (l'^n(fc)l'ffnK + iWn{k)?9lK) , (3-8) 

nfeo 
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where the subscripts s and t stand for 'singlet' (j = 
0, m = 0) and 'triplet' (j = l,m = ±1,0 ) respec- 
tively; the factor 3 in front of g^^^ accounts for the triplet 
degeneracy. For a one-dimensional system, the exciton 
wave function near band buttom can be approximated 

byl(/)J*(fc) = i2r^''/7ry/y{l+[ik±n/2y/f}, where 
Tq* are the radii of the excitons. Note that the conduc- 
tion band bottoms are at T'''/2- The ground state makes 
major contribution since all the other levels are far off 
resonance and barely occupied. According to Abe's sin- 
gle configuration interaction (SCI) calculation,E3 the ra- 
tio Uq/Vq determines whether the singlet exciton level 
( ^Bu) or the triplet exciton level ( is lower in en- 

ergy. They found that, for Vb = to and Uq > 1.39Vb, 
is lower than and vice versa. All the even-parity Ag 

states are lying at higher energies. However, it is found 
in some finite-chain calculations that-the lowest excita- 
tion is actually an even parity state.Ej This would have 
a significant effect on the efficiency of luminescence since 
the optically excited electrons at ^Bu may relax to the 
Ag state first, then release their energy via non-radiative 
channels. Nontheless, it has been shown that the relax- 
ation rate from the 1 ^B^ exciton to the 2 ^Ag exciton 
is much smaller than the relaxation rate to the ground 
state.E^ Furthermore, from the point of view of phase 
space filling, both 1 ^i?„ and 2 ^Ag excitons contribute 
almost equally to the reduction of phase space. Both 
factors seem to diminish the effect of this even-parity 
state. In fact, for the optical nonlinearity being studied, 
one study shows that, for both SCI and finite-chain cal- 
culations, the optical nonlinearity is determined almost 
entirely by the odd-parity ^B^^ exciton, a dominant Ag 
exciton above the ^Bu level (not considered here), and 
the threshold of the conduction bandH3 Therefore, as far 
as optical nonlinearity is concerned, we will neglect the 
influence of this even-parity state in this paper. 

The thermal distribution of excitons is given by g^'^ = 
(exp{/3[e^-*(i^) - ^,,«]} - where ef^\K) = e^'*(0) + 

/2mex are the energies of excitons, the exciton effec- 
tive mass rriex is fi^ /o^tQ, and ^s,t are the quasi-chemical 
potentials. Immediately after the optical pumping, there 
are only singlet excitons because of the selection rule. 
Part of these excitons then fall down to the triplet level 
via spin-orbital interaction. Their populations are con- 
trolled by the quasi-chemical potentials in our calcula- 
tion. 

After summing over all the electrons in the conduction 
band, we have 

fci ^ 



citon wa ve functions. By taking the expectation value of 
Eq. (3^); we have 



2MNc 

nn'jm kK 

± \^ t 

^ njmK 



5Z \^njmK^njmK 
njmK 



(3.10) 



(3.9) 



where we have used the completeness relation for the ex- 



This identity is used to determine the values of the chem- 
ical potentials /is and fit, once the total population of, 
and the relative populations between singlet and triplet 
excitons are given. 



IV. ABSORPTION SPECTRA: NUMERICAL 
RESULT 

Most of this section is devoted to the calculation of 
the resonant nonlinear optical spectra of PDA for reasons 
stated at the introduction. At the end of this section we 
will comment briefly on the calculation for PPV. 

PDA has four carbon atoms per unit cell, and conse- 
quently four bands in the tight-binding approximation. 
Since our focus is on the exciton state within the band 
gap, which is composed of the electron and hole from the 
middle two bands (conduction and valence band), the 
outer two bands can be safely neglected. We choose the 
dimerization-induced band gap, ASt, to be the unit of en- 
ergy. The value of to, which determines the total band 
width oLconduction band and valence band, is chosen to. 
be 1.250 The average bond length of PDA is 1.35 A.eI 
The Coulomb interaction parameter Vq = e^/eao can 
be determined from the intra-chain dielectric constant 
e. According to the Kramers-Kronig analysis of the re- 
flectivity and-witb permittivity measurements, this value 
is close to 3.E2I We choose e = 3.5, such that Vq = 2.84 
eV, to fit the calculated exciton binding energy with the 
value observed in experiments. In most |Qf the following 
calculations the on-site energy Uq = 2Voc3 We have also 
done several calculations using different choices of Uq and 
Vo values, (see Fig. 4) In all of the following calculations 
the number of sites is 400 and the damping 7 is 0.02. We 
have done calculations on a larger system with 800 sites 
and confirmed that the finite size effect is unimportant. 

Figure 1 displays the electric susceptibility x(w), in 
room temperature, for various concentrations of the sin- 
glet exciton gas. The percentages of electrons excited 
from the valence band by the pump wave are indicated 
in the legend. The absorption peak in Im x(w) can be 
clearly identified at huj = 1.39 in the absence of pump- 
ing, (the solid line) The rcnormalized band gap is de- 
termined by the first maximum of Im x(^) beyond the 
exciton peak, which is barely observable at tko = 1.71. 
Notice that the band gap renormalization due to the 
Coulomb interaction is quite large. With different ex- 
citon populations, the magnitude of renormalization also 
changes, (see the discussion after Eq. (2/7)) The exact 
value of the conduction band edge is difficult to measure 
experimentally because most of the oscillator strength is 
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'concentrated' on the exciton peak. On the other hand, 
the singlet exciton absorption peak at 1.97 eV is one of 



jAficurately and is 
BEa Therefore, it 



the few values that can be determinec 
generally agreed upon by researchers .L 
is used to set the overall energy scale and that means 
1.39 X iSt = 1.97 eV, or iSt = 1.42 eV. Consequently, 
the position of conduction band edge is at 2.43 eV in our 
calculation and the binding energy of the singlet exciton 
is approximately 0.46 eV. This_Ealls within the range 0.4 
- 0.5 for the values reported.Elllj 



concentration increases. The reduction is approximately 
proportional to the number of excited electrons- This is 
consistent with the picture of PSF, which givesEi 



Sf/f 



sat 
ex : 



(4.1) 



where / is the oscillator strength, Uex is the exciton den- 
sity per unit length, and n*"* is the saturation density, at 
which exciton wave functions begin to overlap in space. 
By using Eq. (4.1), we can estimate the exciton radius 
and the result is about five unit iseUs, which is again con- 
sistent with earlier calculations.tj't3 
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Fig 1: Real (a) and imaginary (b) parts of x(w) for dif- 
ferent exciton concentrations, labeled by the percentages 
of the occupation of conduction band (0, 0.5%, 1.0%, and 
1.5%) at room temperature. Probe energy is in units of 
the bare band gap, A5t. The 'vvcc' curve comes from the 
calculation without the V^^^^ and V^l^^^ terms and in the 
absence of excitons. The arrow at the bottom right of 
Fig. 1(b) indicates the position of the band edge for the 
solid line. 

In Fig. 1, the absorption signal is reduced as exciton 



Fig 2: The effect of different exciton concentrations on 
(a) absorption cocfhcient a and (b) the change of refrac- 
tion index 5n. a is in units of 2.0 x 10^ cm^^. The energy 
is in units of ASt. 

Figure 2 shows the absorption coefficient a and the 
change of refraction index 5n near the exciton peak. The 
values of the exciton concentration are the same as in 
Fig. 1. The chain-to-chainj^istance for PDA in the crys- 
talline phase is about lOAJj and the inter-site distance is 
I.4A. Therefore, 1 % concentration has 7.1 x 10^^ cm^"^ 
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conduction electrons. The pump intensities Ip required 
for this electron concentration be obtained from 

Ip = TiexfiLOp/ oitT (reflection from the sample is ignored), 
where TiLjp — 1.97 eV, the peak absorption at (includ- 
ing 7r-eleC|t|i'ons and the background) is approximately 
10^ cni^^,Q and the recombination time r is 2 ps.El Con- 
sequently, to excite 7.1 x 10^^ cm~'^ electrons requires a 
pump wave with intensity Ip = 1.14 x 10^ W/cm^. The 
optical Kerr coefhcient n2, which measures the change 
of refraction index due to pumping, is given by n2 = 
\5n\/5Ip. For a pulse with Ip = 1.14 x lO'^ W/cm^, we 
have \5n\ — 0.196 at resonance (see Fig. 2b) and therefore 
712 = 1.7 X 10^^ cm^/W. This is four orders of magnitude 
larger than the nonresonant value, ancLis close to Greene 
et al's observation 3.0 x 10"*^ cm^/W-H 
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Fig 3: The effect of temperatures and singlet /triplet 
populations on Im x(w). The percentage of conduction 
electrons is 2.5% for all of the curves. At room temper- 
ature, it makes small difference whether the excitons are 
singlet or triplet. The distinction is more apparent at 10 
K. 

All of the above calculations are for singlet excitons at 
room temperature. Using the same formalism it is quite 
easy to investigate the influence of exciton species and 
temperature on the resonant nonlinearity. We consider a 
two- level model where only the singlet ( ^Bu) and triplet 
( '^Bu) excitons are consideredj-The exciton radii are cho- 
sen to be 6a (rg) and 4.5a (?'o)ilj In Figure 3, we show the 
extreme cases where the populations are either all-singlet 
or all-triplet. This difference traces back to the different 
distributions of the relative part of the electron-hole pair 
wave functions (/)Q'*(fc). It can be seen that the difference 
between singlet and triplet curves are more significant at 
low temperature (lOK). Such a temperature effect has 
not been studied experimentally, however. 

In Fig. 4, we show how different choices of the strength 
of the electron-electron interaction may affect the absorp- 
tion spectra. It can be seen that the magnitude of Vq, the 



long-range interaction, has significant effect on the band 
gap renormalization and the binding energy. For exam- 
ple, the band edges for Uq — i and Vq— 1.5, 2, 2.5 are 
at 1.54, 1.71, and 1.89, respectively (see arrows in the 
figure). The widths of band gap vary roughly linearly 
with Vq. On the other hand, the short-ranged on-site 
energy Uq has little influence on the bandgap. Further- 
more, contrary to the effect of Vq, a larger Uq leads to a 
smaller binding energy. This adverse effect can be traced 



back to the exchange terms ^x'^x^^s^i ^"i' (^)- 



dependence of the position of the exciton level on Uq 
Vq resembles closely to the calculations by Abe et 
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Fig 4: The influence of Columbic parameters, Uq and 
Vq, on Im xi'-^) in the absence of pumping. The value 
of {Ua — Vq) is indicated for each curve. The arrows 
(from left to right) indicate the band edges for the curves 
(4-1.5), (4-2.0), and (4-2.5) respectively. Their po- 
sitions shift linearly with Vq. On the contrary, the band 
edge is not affected by varying Uq (cannot be seen easily 
from the flgure). The energy is in units of ASt. 

This work is based on the simplifled model of an infi- 
nite and rigid chain and seems to be applicable to other 
conjugated polymers, wjth exciton levels as well, such as 
PPV and polysilane.Bo At the end of this section, we 
present a calculation for PPV. Its optical nonlinearity 
(for a perfect crystalline sample) is found to be of the 
same order as PDA's. The following parameters obtained 
fromptiie spin density profile of nonlinear excitation are 
usedfea to = 2.02 eV, Uo = 2.5to, K, = l-3io. The dis- 
placement of ions due to the double bond alternation 
is about ±0.055A. This gives ASt — IJ5 eV because the 
lattice stiffness for PPV is 5.23 eV/AH3 From the absorp- 
tion spectrum being calculated, we estimate the exciton 
binding energy tOr-be 0.6 eV, wiiile earlier calculations 
range from 0.4 eVE2l to 0.8 eV.Ea The exciton radius in 
PPV is about 50 A, which spans 8 unii-cells. The inter- 
chain distance for PPV is about 4 Af3 the absorption 
coefficient for PPV at 400 nm is 2.3 x 10^/cm at room 
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temperature, and the total exciton lifetime for PPV film 
is 0.32 ns.EZI From these data and the change of refrac- 
tion index being calculated, we find that the optical Kerr 
coefHcient for PPV to be 8 x IQ-^crnVW. This is of the 
same order as the value for PDA. 

However, in actucal practice, PPV is rarely used in the 
study of optical nonlinearity, whether it is resonant or off- 
resonant. This may be related to engineering or chemical 
problems in growing high quality single crystals. These 
aspects are beyond the scope of this paper. Even if a high 
quality single crystal for PPV can be obtained, the inter- 
chain distance will be much smaller in PPV (3 to 4 A, 
comparing to 10 A for PDA). Under this circumstance, 
the effect of inter-chain coupling probably will invalidate 
our result presented here. Another complication for PPV 
is that, for the intensity of 10^ W/cm^, we probably have 
entered thejregime where exciton-exciton annihilation is 
significant .EZl Also, because of the high concentration the 
possibility of excited-state absorption increases. Thus, 
actual pumping efficiency should be lower than that re- 
ported here. All of these reasons will more or less make 
our calculation futile for PPV. 

On the other hand, exciton-exciton annihilation may 
not play a significant role in PDA below pump inten- 
sity 10* W/cm^. A recent paper by Schmid showed that 
the susceptibility of PDA does not deviate fromjA pure 
X^'^) behavior until the peak intensity 10* W/cm^.E3 This 
seems to indicate that, neither exciton-exciton annihila- 
tion, nor excite-state absorption, is appreciable within 
the range of our consideration. One possible reason for 
the higher threshhold of exciton-exciton annihilation is 
that the exciton radius in PDA (about 10 A) is much 
smaller than that in PPV (about 50 A) and has higher 
saturation density. This explains to some extent why the 
same formalism works so well for PDA, but not for PPV. 
This is also supported by the fact that simple estimates 
on the exciton density without such annihilation effect 
has been quite consistent with the experimentsQIa. 



V. CONCLUSION 

Resonant optical nonlinearity for conjugated polymers 
can be understood using the simple picture of PSF, in 
which the probe signal is reduced because the phase space 
for final states has been occupied by the excitons. This 
model provides only an order of magnitude estimate and 
fails to produce more details such as the probe response 
over the whole spectral range (eg. the band edge ab- 
sorption), the position of the exciton peak, the effect 
of temperature and Coulomb interaction... etc. The 
present study is able to access the effects of various micro- 
scopic parameters by using the EOM method. We pro- 
duced the electric susceptibility xi^) that contains the 
information about the positions and oscillator strength 
of exciton level and conduction band edge. By varying 
the exciton populations, we can observe the change of 



the resonant oscillator strengt h an d the trend is consis- 
tent with the PSF model (Eq.(^). We also calculated 
the optical Kerr coefficient n2 and the value obtained 
1.7 X 10~* cm^/W agrees well with observations. 

It has to be reminded that several complications in a 
real polymer system have been left out to simplify the dis- 
cussion. We have used a rigid and infinite polymer chain 
while in actual experiments it is finite and maybe flexible. 
The phonon degrees of freedom will contribute to extra 
features,in the absorption curve such as the phonon side- 
bandsEj We have also used the quasi-equilibrium condi- 
tion for the exciton gas. In future research, the dynamical 
evolution of exciton density can be included by coupling 
the equation of motion to the rate equation of the elec- 
tron population. Finally, the present theory has to be 
modified at high exciton density when exciton-exciton 
interaction plays a more important role and the RPA is 
no longer valid. This may lead to exciton-exciton anni- 
hilation, formation of bi-excitons, and even the existetice 
of a gain threshold beyond which lasing can happen.li2l 
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APPENDIX A: 

Based on the symmetry of inversion, we can sh ow t hat 
the terms that are quadratic in Uc^vka in Eqs. ( ^.6[ ) do 
not affect the dynamics of the total polarization. First, 
the inversion properties of the Bloch states are, 

*c/c(-7-) = -«'c-fe(r), 

^vk{-r) = ■i'v^kir). (Al) 

This leads to the following identities: 

lJ-vc{-k) = Hvc{k), (A2) 

and 

"'^AiAaAaAi ^ "^^AfAj AsAi^'^i ""As ^TAg '''Ai : 

^AiaIa3A4 = "^AfA2A3A4'^Ai7rA2 7rA3 7rA4, (A3) 

where tTc = — 1, 7r„ = 1. It is clear that '^x'^x^x^^x^ changes 
sign when the number of c indices does not equal to the 
number of v indices. 

Second, combining the exciton wave function, 
</)^^*(fc) = (2r^'77r)VV{l + [(fc±7r/2)r^'f }, and 
the Bose-Einstein distribution function, g'^^j^ — 
{exp{(3[ef;\K) - Hs,t]}- 1)"^ (see Sec. HI B), and using 
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the fact that ^'fe(r) = ^'/c±7r(r), it is not difficult to see 
from Eq. (3.8) that 

Because of the symmetries in Eqs. (A2), (A3), and 
(A4), the equation for p^k has the same form as the equa- 
tion for pk (see Eq. ( ), with k and k' being replaced by 
— fc and —k') except that the signs of the terms quadratic 
in n-c^vkcr are changed. Consequently, they do not con- 
tribute to the total polarization (P), in which pk and 
P-k appear through the combination of iivc{k){pk +P-k) 
only. Notice that the conclusion may not be valid if the 
electron population is not in thermal equilibrium. 
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